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The statistical mechanics of quarks and gluons are investigated within the context of the canon- 
» I ■ ical ensemble. Recursive techniques are developed which enforce the exact conservation of baryon 

number, total isospin, electric charge, strangeness, and color. Bose and Fermi-Dirac statistics are 
also accounted for to all orders. The energy, entropy and particle number densities are shown to be 
significantly reduced for volumes less than 5 fm 3 . 
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> ' I. INTRODUCTION 

m ' 

The study of the quark gluon plasma in the laboratory through hadronic collisions represents one of the largest 
■ current initiatives in nuclear physics. Although these volumes are manifestly finite, modeling of the collisions seldom 
addresses the impact of the finiteness as most treatments are based on the the grand canonical ensemble where the 
Cf") . conservation constraints of baryon number, strangeness, isospin and color are ignored. Indeed, one expects that these 
constraints would become meaningless for larger volumes. The effective volume for a pp collision might be only a few 
cubic fm, while the volume of a RHIC collision might be a few thousand cubic fm. But even in a RHIC collision, the 
duration of the quark gluon phase might be only a few fm/c, meaning that the effective volume for charge conservation 
may be a few dozen cubic fm. 

Recently, the effects of local charge conservation have gained more attention due to the relationship of charge 
conservation to charge fluctuations and charge balance functions 0, 0j|E0,IEIE0I1)@1- These observables are 
intimately related to the delay of hadronization expected should a quark-gluon plasma be created. Recent papers 
have focused on the importance of conserving electric charge, strangeness and baryon number, as well as the dynamics 
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of charge equilibration [ljj . The statistical mechanics of a finite- volume hadron gas have been studied in the context 
j_j ■ of pion flavor distributions [TTj . 

The approaches for solving the canonical ensemble can be divided into two classes. In the first set of approaches, one 
can sum over all partition functions in the grand canonical ensemble assuming the chemical potentials are imaginary, 
Z ~ Tre~' 3ff+l ' 3A ' < 2. By integrat ing over all (i, the phase factor becomes a delta function which restricts the phase s pac e 
to states with Q = |l%llall4| . This approach can even be applied to non-Abelian symmetries such as SUfN) |l5lll6| . 
By associating N— 1 charges, Q\-Qn-i, with the N — 1-fold Cartan subgroup, one can restrict the ensemble to states 
with specific charges Q, which can then be associated with the SU(3) multiplet labels, (p,q), by a transformation. 
These techniques work well for simple systems, e.g. massless partons, where the integral over complex phases can be 
be performed analytically. Multiple charges can be readily incorporated, including conservation of overall momentum. 

The second class of approaches for the canonical ensemble centers about recursion relations which were proposed 
by Chase and Mekiian|l7l for the study of nuclear fragmentation. These methods have been extended to include 
quantum statistics |19| and non- additive charges , such as angular momentum and isospin. In addition to partition 
functions, these methods can also generate multiplicity distributions [1 ll |2C| . The advantage of these methods is that, 
given the one-particle partition functions, they provide exact answers. Since these methods involve sums, rather 
than integration over phases, they provide robust answers for arbitrary energy levels. Furthermore, symmetrization 
is included to all orders. If the number of particles to be considered remains below 100, numerical calculations tend 
to take only seconds or minutes at most. 

We review recursive techniques for canonical ensembles and present extensions to incorporate conservation of color 
in the next section. In the subsequent section a simple example of a non-interacting quark-gluon gas is explored. 
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Canonical partition functions are generated which enforce conservation of baryon number, strangeness and charge, 
while requiring the system to be in an isosinglet and in a color singlet. Bose and Fermi statistics are included to all 
order. Energy, entropy, and particle number densities are all shown to be significantly suppressed for small volumes, 
less than 5 fm 3 . Arguments are made that these effects are unlikely to be important in a relativistic heavy ion collision, 
but might be quite important for pp collisions. 



II. RECURSIVE TECHNIQUES FOR QUANTUM COLORED PARTICLES 

For additive charges, canonical partition functions can be calculated in a straight-forward manner. Neglecting 
symmetrization, conservation of a set of charges Q can be enforced through the recursion relation, 

^d T )=J2 aj ^ z A- ah ,^(n (i) 

Here, k refers to all particles species of charge q k whose single-particle partition functions are: 

^ = 5> e - £ */ T , (2) 

i 

where i denotes the individual energy levels of energy Ei and degeneracy di . The number a k must be positive for every 
species. For low energy nuclear physics applications, a k could be the baryon number and A would represent one of the 
conserved charges. For a high-energy system where ant i- particles must be included and there are no positive-definite 
conserved charges, A can refer to the number of hadrons. This number is positive-definite, but as the number of 
hadrons is not conserved, one must consider all A to obtain the partition function, Z = Za- 
Symmetrization can be incorporated by considering all permutations of identical particles |l9t l2l| , 



Z a,q - Al 



i E ^■■■jA\e- H l T \V{j l ...j A )){-ir^\ (3) 



where a sum over all Al permutations of the single particle levels accounts for the symmetrization and Np(V) counts 
the number of pair- wise interchanges of fermions involved in the permutation V. 
By considering the cyclic permutation, 



c k,i {T) _ J- {juh . ..jA e - H / T \j 2 , ■ ■ 

h—h 



E< 



(4) 



one can derive a modified recursion relation, 

~A 



Z A^) = EE^(±1^^P0. (5) 



k i 



By summing over all possible values of £, symmetrization is included to all orders. 

Before addressing SU(3) symmetries, we review the techniques used for SU(2) as described in ^lj. Extending 
the techniques described above to incorporate an SU(2) symmetry such as isospin or angular momentum can be 
accomplished in two ways. First, one can consider the projection M as an additive charge and solve for Zm(T) with 
the method outlined above. Given the (2/+ 1) degeneracy of all multiplets of size /, the partition function for states 
with fixed I is easily obtained from the partition function for fixed M, Zj = Zm — Zm+i- The second technique 
requires one to understand the decomposition of the cycle diagrams into SU(2) multiplets. In terms of the projection 
operator, which projects states with total isospin I and projection M, one can define the isospin decomposition 
of the cycle diagram. 

Cf{T) = E (h,h ■ ■■3i\e- H/T O i m \h, ■ ■ ■ 

ji---jt (6) 

=u K {T/l) X f, 
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where xT is the operator Cf.' 1 evaluated with the particles being combined to one quantum state of zero energy. The 
subscript k refers to a specific multiplet, not the individual states. For instance, k might refer to pions and ui K {T) 
would be the partition function for a given species of pion. The details about the energy levels are absorbed in uj k 
while x is independent of the energy levels, but does depend on the isospin of the particular species, i K . 
Once Xi' 1 is known, one can calculate the partition function for fixed / recursively. 

K I 

One could have added indices M, M' and m onto Za.i, Za'j 1 and xt' respectively. However, as each of these 
functions involve a trace over all states, and the Hamiltonian has no dependence on the isospin projection, none of 
these quantities has any dependence on the projection and the indices can be suppressed. The function Af(I',i,I) 
counts the number of multiplets of size I that one obtains by coupling /' with i. For SU(2), these obey a simple form. 



Af(l',i;I)= ]T (I,M\I',M',L 



(8) 

i, \r -i\<i <r +i 

0, otherwise 

Before performing the recursive calculations described in Eq. 0, one needs a simple expression for xt' ■ This can be 
accomplished by realizing that by summing over i one obtains the projection of states with fixed m. 

zZ*h'™ = C?" 1 ' & " ' 31 1 °rn \h ■■■U,3l)- (9) 

i jl—jt 

Since the states \ji ■ ■ ■ ji) are eigenstates of O m , the bra and ket must be identical, i.e., all particles must be identical 
and have the same isospin projection. 

sp i,i = f 1, m = 0,±l,±2£,---±B K 

/_^XK,m I otherwise ^ ' 

i 

Since Xi'm ^ s independent of m when m is a member of the isomultiplet i, the sum over i in Eq. (|10|l can be limited 
to i > m. It is then straight-forward to see that 

1, i = 0,£,2i,---U K 
-1, i = I - 1, 21 - 1, • • • li K - 1 (11) 
0, otherwise 

Here, the index m is suppressed since x is independent of m as long as m is a member of the isomultiplet i. For the 
i K = 1 case, e.g., pions or p mesons, x\ — 1 f° r i = £ an d i = 0, and x\ — ~ 1 f° r i = — 1. 

Restricting SU(3) color can be accomplished with the same procedure. Following the same steps as the SU(2) case, 
one can derive analogous expressions, 

Z a,q.(p, Q) (T) = £ °«"«CT/l) (±l)^iz A _ fa|tt ^ t(p<tq<) (r)x^((^, Q'), (p, g); (P, Q)). (12) 

K,l,{P',Q),(p,q) 

Here, (P,Q) denotes an SU(3) multiplet, e.g., the gluon color octet is represented by (P — l,Q = 1). For an 
explanation fo the notation see [22^. Again, the challenge in making Eq. (|12|) tractable is in finding expressions for 
X^p ? ) and N{{P' ', Q'), (p, q); (P, Q)). Coupling SU(3) multiplets, i.e., finding expressions for Af, follows rules based on 
manipulating Young tableaus |22j |. For the calculations in this paper, these rules were programmed numerically. 

Finding an expression for x*p 9 ) requires performing a color decomposition of the cycle diagrams and can be done 
analogously as was shown above for SU(2). First, one must project states of a given /i, where /i represents the 
eigenvalues of the projection operators in SU(3), e.g., hyper-charge and I3 for SU(3) flavor. As with the SU(2) 
example above, this projection is realized by summing x^L q \ ^ over all (p, q) multiplets which include the projection 

X( P e ,q),v = zZ ^i>h ■ --kP^h ■ ■■ii.ii)- (13) 

(p,q) h—h 
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Again, the fi dependence in x represents the known degeneracy of fj, within a given (p, q) multiplet. 

At this point, we proceed by considering gluons as an example, and drop the index k. For £ gluons in a given 
quantum state the trace of (ji,j2 ■ ■ ■ ji\0^\j2 ■ ■ ■ will be zero unless \i corresponds to £ identical gluons. The 

upper-left panel of Fig. ^ displays the projections \i that result from this trace applied for five gluons. The upper- 
right and lower-left panels display the (£,£) and (I — 2,1 + 1) multiplets. By careful inspection, one can see that 
the states along the diagonal, as shown in the lower-right panel of Fig. ^ are represented by the combination, 
{£, £) - {£ + 1, £ - 2) - {£ - 2, £ + 1) + (£ - 1, £ - 1). By subtracting the same object evaluated with £' = I - 1 and 
adding two color singlets, one can see that 

r 1, (p,q) = (£,£) or (£-3, £) or (£,£-3) 
v giuonV_l -1, (p,q) = (£-2,£+l)or(£ + l,£-2)or(£-2,£-2) 

X (p<i) \ 2, (p = 0,g = 0) 1 > 

y 0, otherwise 

Inserting this expression into Eq. (|12|l allows one to account for Bose effects in the gluonic partition function. 

Quarks are represented by the multiplet (p — l,q = 0) while anti-quarks are represented by the multiplet (0, 1). 
Following the same ideas as were illustrated for the gluons, one can show 



1, (p, = 0) or (f- 3,0) 
T, (p,q) = (£ 
0, otherwise 



X?;T' e ={ -1, (p,q) = (£-2,l) (15) 



By combining the recursion relations for additive charges, which are used to account for baryon and strangeness 
conservation, with the recursive techniques for SU(2), which are used to enforce conservation of isospin, and the 
methods for SU(3) which account for color, one can construct a recursive prescription which account for all the 
charges in a parton gas. One might also choose to write partition functions for subsystems, e.g., the strange quarks 
or the gluons, then convolute the partition functions together to find the partition function of the combined system. 
For instance, partitions for the two subsystems a and b can be combined to find the combined partition function, 

z (p,q)( t )= E ^ , fl .)( T ) z L»)( T )- A/ '((p«'««)' (p»'*)"> ( 16 ) 

III. RESULTS AND CONCLUSIONS 

The methods described in the last section were applied to the example of a parton gas in a volume V at a 
temperature, T = 250 MeV. The gluons, up-quarks and down-quarks were assumed to be massless, while the strange 
quark was assumed to have mass of 150 MeV. The single-particle partition functions were calculated by integrating 
over the momenta, 



2V 
mV 



(17) 



(0mK o (0m) +2Ki(j3m)) 



The two helicities were accounted for by the factor of two preceding the expression. One could easily account for shell 
effects by replacing the integral over momentum with a discrete sum. Although discrete states are more consistent 
given the Bose and Fermi effects described above, the continuous form for oj is used here so that charge conservation 
effects can be viewed separately. 

The calculations were performed according to the following prescription: 

1. The gluon partition function, Z|^ u °", was calculated recursively by first calculating for all numbers of gluons A, 
then summing over A. 

2. The partition function for strange quarks, Z S A ^ ^ , was generated. The partition function for anti-strange quarks, 
Z S A ^ s, was then generated by switching (jp,q) with (q,p). The partition function of the strange/anti-strange 
quark system was found by summing over all partition functions with equal numbers of s and s quarks. 

z 'S, q )( T ) = E ^,(p^ 9a )( T )^l(p 6 ,, 6 )( T )^((Pa,ga),b b ,g 6 );b,9)). (18) 

A,(Pa,q q )(pb,qa) 
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3. The partition function for up and down quarks were calculated separately. In terms of the isospin projection, 
m, and the net number of up/down quarks, A, the partition function for the up and down quarks is: 

Ztm, M (T)= J2 Z ™ + A/2Xp , 9o )( T )^^n + A/2,( P6 , 96 )(W(p a , gQ ),(p b , (? ,);(p, (7 )). (19) 

4. The partition function for ud quarks is calculated by interchanging p and q. In the same manner as the ss 
quarks were convoluted, the ud and ud partition functions were convoluted to find the partition function for u, 
d, u and J quarks with the constraint of zero baryon number, Z^^Z s, where M is the isospin projection. 

5. The partition function for the udud system constrained to an isosinglet is calculated by taking the difference of 
the M = and M = 1 partition functions. 

6. Using Eq. i|16|) . the partition function for the udsuds system were generated. This was then convoluted with 
the partition function for gluons to find the partition function for the entire system. 

Energy densities can be calculated with the well-known formula, (E) = {—d/d[3) ln(Z). This requires performing 
the calculation at two adjacent temperatures, doubling the CPU time of calculating the partition function. Thus, 
an alternative approach was developed where recursion relations were derived for {d/d0)Z. Recursion relations 
were also generated for the trace of iVe _/3ff , where N could be one of several number operators, iVgi uons , N s , N u 
or Nd- The recursion relations for these functions were calculated simultaneously with the recursion relations for 
the partition functions with little penalty in CPU time, because all calculations involve the same calculations of 
M(( P ,q),(P',Q');(P,Q)). 

Figurc|2]shows the energy and entropy densities as a function of the volume. Both the entropy, S = \n(Z)+f3(E), and 
the energy are lowered by the reduction in the available states due to the conservation constraints. The grand canonical 
limit, represented by the dashed lines, is approached at large volumes. As shown in |l 5l . restricting the matter to a 
color singlet reduces the partition function by a factor which scales as V~ . The other charge conservation constraints 
result in additional reductions to the quark sectors. A reduction of V~ A to the partition functions corresponds to a 
reduction of — 41n(V) to the entropy and energy densities. Since the bulk contributions scale as V the logarithmic 
conservation penalties become irrelevant at large volumes. The chemical decomposition of the various species behave 
similarly as shown in Fig. [31 The color penalty for gluons is more severe than the color penalties for quarks since 
gluons effectively have larger color charges and are less likely to couple to a singlet. 

From Fig.s|21and|21 it is clear that conservation rules are important for volumes of 5 fm 3 or less. It is therefore 
important to take such considerations into account in pp or pA collisions. Heavy ion collisions at RHIC occupy 
thousands of cubic fm. However, the local nature of charge conservation results in an effective volume which is 
determined by the initial conditions, and the diffusion of the various charges. Given that the color of a projectile 
nucleon is spread over several units of rapidity by the initial stopping process, and given the spread of charge over a 
few fm in the transverse direction, it would be difficult to make a case that the effective volumes are less than a few 
dozen cubic fm. It should be emphasized that the determining factor for the importance of conservation constraints 
is the number of partons in the effective volume. For massless partons, the density of partons is expected to exceed 
10 fm~ 3 in a during the first one or two fm/c of a RHIC collision. If the degrees of freedom are restricted by another 
means, e.g., a large effective mass for partons, the relative penalty for color conservation would increase. 
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X12344321X XXXX1 1 XXXX 
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FIG. 1: Upper-left: A diagrammatic representation of the cyclic trace of five gluons in a single quantum level. Upper-right 
and lower-left: A diagrammatic representation of the (5,5) and (6,3) states. The integers refer to the degeneracy of states 
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for £ — 4, and adding two color singlets, one obtains the desired combination of states in the upper-left diagram. 
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FIG. 2: The energy and entropy densities are shown for a parton system at a temperature of 250 MeV as a function of 
the volume. For small volumes, conservation of baryon number, isospin, strangeness and color restrict the phase space and 
significantly lower the entropy and energy density. For larger volumes, the results approach the grand canonical limit (dashed 
lines). 
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FIG. 3: Densities of gluons (circles), s and s quarks (triangles), and u,d,u and d quarks (squares) are suppressed for small 
volumes due to charge and color conservation. The grand canonical limits are represented by dashed lines. 



